arXiv:hep-th/0503047vl 4 Mar 2005 


Universal Crossover in Perturbation Theory with a 
Large Field Cutoff 


L. Li| and Y. Meuricej |§ 

t Department of Physics and Astronomy 
The University of Iowa 
Iowa City, Iowa 52242 
USA 

I Also at the Obermann Center for Advanced Study, University of Iowa 


Abstract. For models, the introduction of a large field cutoff improves 
significantly the accuracy that can be reached with perturbative series but the 
calculation of the modified coefficients remains a challenging problem. We show 
that this problem can be solved numerically, and in the limits of large and small 
field cutoffs, for the ground state energy of the anharmonic oscillator. For the two 
lowest orders, the approximate formulas obtained in the large field cutoff limit 
extend unexpectedly far in the low field cutoff region. For the higher orders, the 
transition between the small field cutoff regime and the large field cutoff regime 
can be described in terms of a universal function. 
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1. Introduction 

Perturbative methods and Feynman diagrams have played an important role in the 
development of quantum field theory and its applications. However, perturbative series 
usually have a zero radius of convergence For scalar models with A(/)^ interactions, 
the coefficients of perturbative series grow factorially. For any fixed, strictly positive, 
value of A, there exists an order beyond which adding higher order terms diminishes 
the accuracy. This feature will restrict our ability to perform high precision tests of 
the standard model (for instance, g — 2 oi leptons and the hadronic width of the 
during the next decades. 

The large order behavior of the series is dominated by large field configurations 
which have little effect on low energy observables. Introducing a large field cutoff 
HEl in the path integral formulation of scalar field theory, dramatically improves the 
large order behavior of the perturbative series. In two non-trivial examples [3, this 
procedure yields series that apparently have finite radii of convergence and tend to 
values that are exponentially close to the exact ones. This also allows us to define the 
theory for negative or complex values of A, a subject that has raised a lot of interest 
recently H El An important feature of this approach is that for a perturbative 
expansion at a given order in A, it is apparently possible to determine an optimal 
field cut using the strong coupling expansion EICI, bridging the gap between the two 
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expansions. In other words, the modified perturbative methods allows us to take into 
account non-perturbative effects. 

Despite these promising features, calculating the modified coefficients remains a 
challenging technical problem. While developing a new perturbative method (see e.g. 
Ref. [HI ID for the 6 expansion), it is customary to demonstrate the advantages of a 
method with simple integrals and the non-trivial, but well- studiedjTD P, case of the 
anharmonic oscillator. A simple integral has been discussed in Ref. [D- In this Letter, 
we show not only that this program can be completed in the case of the anharmonic 
oscillator, but also that the results show remarkable properties: 


• For the two lowest orders, the approximate formulas obtained in the large field 
cutoff limit extend unexpectedly far in the low field cutoff region. 

• For the higher orders, the transition between the small field cutoff regime and the 
large field cutoff regime can be described in terms of a universal function. 


In the following, we define the model considered, discuss the numerical calculation 
of the modified coefficients, the approximate methods at small and large field cutoff 
and the universal description of the crossover. We use the quantum mechanical 
notations x and uj instead of the field theoretical notations (j) and m. In all numerical 
applications, we use uj = 1. In the path integral formulation, introducing a field cutoff 
means that the paths with |a;(t)| > Xmax for some t are excluded. In the operator 
language, it means that the wave function 'I'(a:) should vanish for \x\ > Xmax- More 
specifically, we will consider the Hamiltonian 


H = \ + V{x) 


with 


V{x) = 




Xx^ 


if 

if 


< 

> a;n 


( 1 ) 


( 2 ) 


Our main interest if the calculation of the modified perturbative series for the ground 
state energy: 


EoiXmax) = ^^E^^\xmax){Xluj'^)’" , 


( 3 ) 




Accurate numerical values for the E^'^ (xmax) can be obtained by combining the 
method of Ref. CH with an expansion in A. The technical aspects will be reported in a 
more detailed article. When Xmax —> oo, we recover the series that has been calculated 
to high order in Ref. CD. In order to allow comparison among the different orders, 
we define the ratios 

Rk{Xmax) = E^^\xmax) / E^\oo) , (4) 


which all tend to 1 in the x^ax oo limit. The numerical values of these ratios are 
shown in Fig. n A striking feature is that the curves for the various orders have 
approximately the same shape. 

The asymptotic behavior of the series can be studied from this data. If Xmax 
is small enough (say below 2 in Fig. then most of the coefficients have Rk ^ 0 
(asymptotic regime) and a linear behavior is observed for ln(|A(^^) when k increases. 
This allows us to infer that the series converges for |A| < with C ~ 65. 
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Figure 1. The numerical values of Rf;{xina,x) for k going from 1 to 10. 


We now discuss an approximate method of calculation of 1 — Rk{xm&x) in the 
large Xmax limit. For k = 0, our expansion parameter will be e = — 1/2) in 

terms of which i?o = 1 + 2e. Assuming a ground state wave function of the form 

^ (x) oc ^ ^ ) (5) 


and expanding at first order in e, we obtain a second order, inhomogeneous, linear 
differential equation for G. Remarkably, this equation can be integrated exactly and 

/*rc py 

G{x) = —2uj / / dze~^^ , (6) 

Jo Jo 

From the boundary condition = 0 we obtain e = —1/G{xmax)- The 

asymptotic expression agrees with the semi-classical estimate of Ref. |3] . We can then 
use our first order result to calculate 

Jo 

The comparison with numerical data is shown in Fig. [3 Unexpectedly, the 
approximation stays accurate in regions where e is by no means small. Integral 
formulas such as Eq. 0 for the excited states will be discussed in a more detailed 
publication. For higher orders, there is solid numerical evidence for a large Xmax 
behavior of the form 


1 - Rk{j 


i DC 

c) CK ^ 


( 8 ) 


We now turn our attention to the opposite limit of small Xmax- The entire 
potential (including the quadratic term) can then be treated as a perturbation. At 
zero-th order, we have the solvable problem of a free particle in a box of size 2xmax- 
Using the conventional techniques of nondegenerate time-independent perturbation 
theory, or numerical method (after a rescaling that allows the use a numerical value 
for the rescaled Xmax), we can calculate the coefficients of the series: 

OO 

m—0 


( 9 ) 
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Figure 2. Numerical values of Ro{xmax) and Ri{xmax) ■ The solid lines 
represent the large Xmax expressions. The broken lines represent lowest orders in 
the small Xmax approximation. 


The lowest order approximations for k — 0 and 1 are shown in Fig. |21 Numerical 
calculations of 6 terms in the series indicate that the series have a finite radius of 
convergence that seem to shrink slowly when the order increases. For fc = 0 {k = 5), 
we apparently have convergence for yCJxmax < 3.3 (^/cJxmax < 2.3). In the case 
fc = 0 or 1, this is sufficient to have a good overlap between the validity range of the 
two limits. However, for higher orders, other methods are necessary to deal with the 
crossover between the two regimes. 

The shape similarities observed in Fig. ^suggest to parametrize Pk in terms of a 
unique, universal, function U, that can be shifted by a fc-dependent quantity that we 
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Figure 3. Rf^{xj^XQ(k)) for fc = 2,... 10 and the universal function U{x) of Eq. 

im . 


denote Xo{k): 

Idk(^X'fnax') — U (^Xmax xo{k)) . (10) 

In order to determine xq, we used the derivative of Rk with respect to Xmax that is bell 
shaped and we parametrized it as . Numerical fits and the generic 

form of the terms of the perturbative series at large Xmax, led us to infer that A = 2lij. 
Adding exponentially small tails on the left (for k large enough) and imposing the 
correct matching on the right, we obtained: 

/ ^max 

da:exp(-2wx^) . (11) 

-OO 

The dependence on the order is encoded in the shift xo(/c). Empirically, they can be 
fitted quite well with 

xo(fc) ~ 0.87+1.13v^ . (12) 

The Rk shifted by Xo(fc) are shown in Fig. 0 One can see that the collapse of the 
various curves is quite good. The function U fits the collapsed data quite well on the 
upper right, but appears to be slightly above the data on the lower left. This can be 
explained by the fact that the function satisfies the duality relation which exchanges 
the small and large field cutoff regions 

U{-x) = l-U{x) , (13) 

which cannot be exact since the approach of 0 when Xmax —^ 0 is power like while in 
the large Xmax limit, the approach of 1 is exponentially small. In these two limits, 
U fails to provide the correct behavior and we need to resort to the approximations 
discussed before. 

We expect similar features in higher dimensional scalar field theory with a UV 
regulator. Generically, a modified perturbative calculation involves a few coefficients 
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with values close to their usual ones, a few coefficients in the crossover region, and 
the rest of the coefficients taking values much smaller than in regular perturbation 
theory. These three regimes are reminiscent of the three regimes encountered when 
calculating renormalization group flows between two fixed points. One should notice 
the similarities between the graphs presented here and those of Refs. [12 El where this 
type of crossover behavior is studied. In the same way, the approximate duality of Eq. 
m reminds of the duality between fixed points of a simplified renormalization group 
equation ini or the exchange between the perturbative and non-perturbative sectors 
of quasi-exactly solvable periodic potentials na We plan to make these analogies 
more specific (or disprove these ideas), by making explicit calculations with Dyson’s 
hierarchical model. Namely, we plan to compare the relative weighs of the three types 
of perturbative terms discussed above with the relative increases in the three regions 
of the renormalization group flows discussed above for the zero-momentum n-point 
functions. 

This research was supported in part by the Department of Energy under Contract 
No. FG02-91ER40664. We thank J. Cook and B. Kessler for valuable discussions. 
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